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Recent theoretical and experimental improvements in the determination of charm and bot- 
tom quark masses are discussed. The final results, me(3GeV) — 986(13)MeV and mf,{my) = 
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986(10)MeV, the presently most precise determinations of these two fundamental Standard 
Model parameters. A critical analysis of the theoretical and experimental uncertainties is pre- 
sented. 
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The past years have witnessed significant improvement in the determination of charm and 
bottom quark masses as a consequence of improvements in experimental techniques as well as 
theoretical calculations. Quark mass determinations can be based on a variety of observables and 
theoretical calculations. The one presently most precise follows an idea advocated by the ITEP 
group more than thirty years ago [|l|], and has gained renewed interest after significant advances 
in higher order perturbative calculations ||^ have been achieved. In particular the four-loop results 
(i.e. the coefficients C„ discussed below) are now available for the Taylor coefficients of the vacuum 
polarization, analytically up to n = 3 and numerically up to n = 10. The method exploits the 
fact that the vacuum polarization function Yl{q^) and its derivatives, evaluated at = 0, can be 
considered short distance quantities with an inverse scale characterized by the distance between 
the reference point = and the location of the threshold q^ = (3 GeV)^ and = (10 GeV)^ 
for charm and bottom, respectively. This idea has been taken up in |^] after the first three-loop 
evaluation of the moments became available ^, |^ and has been further improved in [Q] using 
four-loop results ^ ^ for the lowest moment. An analysis which is based on the most recent 
theoretical | IC , |ri|, |l2|] and experimental progress has been performed in [ 13 1 and will be reviewed 



in the following. 

Let us recall some basic notation and definitions. The vacuum polarization YlQ{q^) induced 
by a heavy quark Q with charge Qq (ignoring in this short note the so-called singlet contributions), 
is an analytic function with poles and a branch cut at and above q^ = Mj^^. Its Taylor coefficients 
C„, defined through 

^Q(^')^QQTLLCnz" (1) 

can be evaluated in pQCD, presently up to order a^. Here z = q^/4nrQ, where mg = is the 

running MS mass at scale jj.. Using a once-subtracted dispersion relation 

I2n^ Jo s[s — q^) 

(with Rq denoting the familiar /?-ratio for the production of heavy quarks with flavour Q), the 
Taylor coefficients can be expressed through moments of Rq. Equating perturbatively calculated 
and experimentally measured moments, 

= I ^^e(^) (3) 

leads to an (?i-dependent) determination of the quark mass 



Significant progress has been made in the perturbative evaluation of the moments since the 
first analysis of the ITEP group. The ff{ccj) contribution (three loops) has been evaluated more 
than 13 years ago [Q, |5|, ^, as far as the terms up to n = 8 are concerned, recently even up to 



« = 30 [ 14, 15]. About ten years later the lowest two moments (n = 0, 1) of the vector correlator 



were evaluated in &{a^), i. e. in four-loop approximation [P, P]. The corresponding two lowest 
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Figure 1: Comparison of resc aled CLEO data for Rb with BABAR data. The black bar on the right 

corresponds to the theory prediction [pO[]. 



moments for the pseudoscalar correlator were obtained in [[T^] in order to derive the charmed 
quark mass from lattice simulations [[l^]. In [|l^, 11] the second and third moments were evaluated 
for vector, axial and pseudoscalar correlators. Combining, finally, these results with information 
about the threshold and high-energy behaviour in the form of a Pade approximation, the full q^- 
dependence of all four correlators was reconstructed and the next moments, from four up to ten, 
were obtained with adequate accuracy [|T2[]. 

Most of the experimental input had already been compiled and exploited in [Q], where it is 
described in more detail. However, until recently the only measurement of the cross section above 
but still close to the B-meson threshold was performed by the CLEO collaboration more than twenty 



years ago [18]. Its large systematic uncertainty was responsible for a sizable fraction of the final 



error on mi,. This measurement has been recently superseded by a measurement of BABAR [ |19| ] 
with a systematic error between 2 and 3%. In [ |T3| ] the radiative corrections were unfolded and used 
to obtain a significantly improved determination of the moments. The final results for mc(3 GeV) 
and m/,(10 GeV) are listed in Table Despite the significant differences in the composition of the 
errors, the results for different values of n are perfectly consistent. For charm the result from n = \ 
has the smallest dependence on the strong coupling and the smallest total error, which we take as 
our final value 



mc(3 GeV) =986(13) MeV, 



(5) 
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n 


mc(3 GeV) 


exp 


a.v 




np 


total 


mfc(lOGeV) 


exp 






total 


1 


986 


9 


9 
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1 


13 


3597 


14 


7 


2 


16 


2 


976 


6 


14 


5 





16 


3610 


10 


12 


3 


16 


3 


978 


5 


15 


7 


2 


17 


3619 


8 


14 


6 


18 


4 


1004 


3 


9 


31 


7 


33 


3631 


6 


15 


20 


26 



Table 1: Results for OTt(3 GeV) and for OTfe(10 GeV) in MeV The errors are from experiment, a^, the 
variation of jX and (for m^) the gluon condensate. 



and consider its consistency with « = 2, 3 and 4 as additional confirmation. Transforming this to the 
scale-invariant mass (m^) pT]], including the four-loop coefficients of the renormalization group 
functions one finds [13] mc{mc) = 1279(13) MeV. Let us recall at this point that a recent study 



[17], combining a lattice simulation for the data for the pseudoscalar correlator with the perturba- 
tive three- and four-loop result [Q |T6[, 11] has led to mc{3 GeV) = 986(10) MeV in remarkable 
agreement with [Q [l3| ]. 

The treatment of the bottom quark case proceeds along similar lines. However, in order to 
suppress the theoretically evaluated input above 11.2 GeV (which corresponds to roughly 60% for 
the lowest, 40% for the second and 26% for the third moment), the result from the second moment 
has been adopted as our final result, 

m;,(10 GeV) = 3610(16) MeV, (6) 

corresponding to mh{mh) = 4163(16)MeV. The explicit dependence of and nij, can be found 



in [|13[]. When considering the ratio of charm and bottom quark masses, part of the and of the jJ, 
dependence cancels 

mc(3GeV) a, -0.1 189 



: 0.2732 - 



mfo(lOGeV) ' " 0.002 
which might be a useful input in ongoing analysis of bottom decays. 



0.0014 ±0.0028, 



(7) 
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Figure 2: Comparison of recent determinations of mc{3 GeV) and miy{mt,). 

In Fig. ^ the results of this analysis are compared to others based on completely different 
methods. The value is well within the range suggested by other determinations. In case of 
our result is somewhat towards the low side, although still consistent with most other results. 
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The results presented in [13] constitute the most precise values for the charm- and bottom- 
quark masses available to date. Nevertheless it is tempting to point to the dominant errors and 
thus identify potential improvements. In the case of the charmed quark the error is dominated 
by the parametric uncertainty in the strong coupling as{Mz) = 0.1185 ±0.002. Experimental and 
theoretical errors are comparable, the former being dominated by the electronic width of the narrow 
resonances. In principle this error could be further reduced by the high luminosity measurements at 
BESS III. A further reduction of the (already tiny) theory error, e. g. through a five-loop calculation 
looks difficult. Further confidence in our result can be obtained from the comparison with the 
forementioned lattice evaluation. 

Improvements in the bottom quark mass determination could originate from the experimental 
input, e. g. through an improved determination of the electronic widths of the narrow T resonances 
or through a second, independent measurement of the R ratio in the region from the T(45') up to 
1 1.2 GeV. As shown in Fig. |l|, there is a slight mismatch between the theory prediction above 1 1.2 
GeV and the data in the region below with their systematic error of less than 3%. 

In this connection it may be useful to collect the most important pieces of evidence supporting 
this remarkably small error. Part of the discussion is applicable to both charm and bottom, part 
is specific to only one of them. In particular for charm, but to some extent also for bottom, the 
/X -dependence of the result increases for the higher moments, starting with n = 4, and dominates 
the total error. We will therefore concentrate on the moments n = \,2, and 3 which were used for 
the mass determination, results for « = 4 will only be mentioned for illustration. 
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Figure 3: mc(3 GeV) for n — 1,2,3 and 4. For each value of n the results from left to right correspond the 
inclusion of terms of order a?, a}, a? and a^. 



Let us start with charm. Right at the beginning it should be emphasized that the primary 
quantity to be determined is the running mass at the scale of 3 GeV, the scale characteristic for 
the production threshold and thus for the process. Furthermore, at this scale the strong coupling 
av(3 GeV) = 0.258 is already sufficiently small such that the higher order terms in the perturbative 
series decrease rapidly. Last not least, for many other observable of interest, like B-meson decays 
into charm, or processes involving virtual charm quarks like B ^ X^J or K ^ nvv, the character- 
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istic scale is of order 3 GeV or higher. Artificially running the mass first down to ^(1 GeV) and 
then back to a higher scale thus leads to an unnecessary inflation of the error. 

The quark mass determination is affected by the theory uncertainty, resulting in particular 
from our ignorance of yet uncalculated higher orders, and by the error in the evaluation of the 
experimental moments. The former has been estimated by evaluating mc(/i) at different renor- 
malization scales between 2 and 4 GeV (changing of course the coefficients C„ appropriately) and 
subsequently evolving mc{li) to mc(3 GeV). The error estimates based on these considerations are 
Usted in Table |I[ 

The stability of the result upon inclusion of higher orders is also evident from Fig. || where the 
results from different values of n are displayed separately in order q'. with / = 0, 1,2 and 3. This 
argument can be made more quantitatively by rewriting eq. in the form 



nir 



4 .^rv ^ 


I (1) 1 (2) 2 , (3) 3 
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OCs - 
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0.618 
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0.299 


\ 0.662 J 




y 0.575 ) 




1^0.396 J 





(8) 



where the entries correspond to the moments with « = 1, 2, 3 and 4. Note, that the coefficients are 
decreasing with increasing order of a^. Estimating the relative error through r™''^a, (3 GeV)'^ leads 
to 1.4 / 2.3 / 2.7 / 2.9 permille and thus to an estimate clearly smaller than the one based on the 
/I -dependence. 

The consistency between the results for different values of n is another piece of evidence 
(Fig. ^ and Table |l|). For the lowest three moments the variation between the maximal and the 
minimal value amounts to 10 MeV only. This, in addition, points to the selfconsistency of our 
data set. Let us illustrate this aspect by a critical discussion of the "continuum contribution", i.e. 
the region above 4.8 GeV, where data points are available at widely separated points only. Instead 
of experimental data the theory prediction for R{s) has been employed for the evaluation of the 
contribution to the moments. If the true contribution from this region would be shifted down 
by, say, 10%, this would move m<^., as derived from « = 1, up by about 20 MeV. However, this 
same shift would lead to a small increase by 3 MeV for « = 2 and leave the results for higher 
n higher practically unchanged. Furthermore, theory predictions and measurements in the region 
from 4.8 GeV up to the bottom-meson threshold, wherever available, are in excellent agreement, as 
shown in Fig. ^ with deviations well within the statistical and systematical error of 3 to 5%. Last 
not least, the result described above is in perfect agreement with the recent lattice determination 
mentioned above. 

Let us now discuss beauty, with evaluated at /X = 10 GeV. Again we first study the sta- 
bility of the perturbative expansion, subsequently the consistency of the experimental input. With 
qs(10 GeV) = 0.180 the higher order corrections decrease even more rapidly. Varying the scale 
jj. between 5 and 15 GeV leads to a shift between 2 and 6 MeV (Table [l|) which is completely 
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Figure 4: R{s) for different energy intervals around the charm threshold region. The solid line corresponds 
to the theoretical prediction. 



negligible. Alternatively we may consider the analogue of eq. (||) with the correction factor 



mb/m\ 



Bom 

b 



/ 0.270 \ 
0.456 
0.546 

\ 0.603 / 



/ 0.206 \ 
0.272 
0.348 

\ 0.410 J 



/ -0.064 \ 

0.048 

0.051 
\ 0.012 / 



a?. 



(9) 



Taking r^'^af for an error estimate leads to a relative error of .28 / .48 / .57 / 0.63 permille for « = 1, 
2, 3 and 4 respectively, which is again smaller than our previous estimate. Let us now move to a 
critical discussion of the experimental input. The contribution from the lowest four T-resonances 



has been taken directly from PDG [ |22| ] with systematic errors of the lowest three added linearly. 



The analysis [ ]13| ] of a recent measurement [ ]19P of Rh in the threshold region up to 11.20 GeV has 
provided results consistent with the earUer analysis but has lead to a significant reduction of the 
error in mt,. 

In comparison with the charm analysis a larger contribution arises from the region where data 
are substituted by the theoretically predicted Rt, with relative contributions of 63, 41, 26 and 17 
percent for « = 1, 2, 3 and 4 respectively. This is particularly valid for the lowest moment. For this 
reason we prefer to use the result from n = 2, alternatively we could have also taken the one from 
n = 3. Let us now collect the arguments in favour of this approach: 

i) For light and charmed quarks the prediction for R based on pQCD works extremely well 
already two to three GeV above threshold. No systematic shift has been observed between theory 
and experiment, in the case of massless quarks, starting from around 2 GeV, and for the cross 
section including charm at and above 5 GeV up to the bottom threshold (Fig. It is thus highly 
unplausible that the same approach should fail for bottom production. 

//) If the true Rb in the continuum above 11.2 GeV would differ from the theory prediction 
by a sizable amount, the results for « = 1 , 2 and 3 would be mutually inconsistent. Specifically, a 
shift of the continuum term by 5% would move mi,, derived from n = 1,2 and 3 by about 64 MeV, 
21 MeV and 9 MeV respectively. 
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To summarize: Charm and bottom quark mass determinations have made significant progress 
during the past years. A further reduction of the theoretical and experimental error seems difficult at 
present. However, independent experimental results on the R ratio would help to further consolidate 
the present situation. The confirmation by a recent lattice analysis with similarly small uncertainty 
gives additional confidence in the result for mc. 

Acknowledgments: I would like to thank K. Chetyrkin, A. Maier, P. Maierhofer, P. Marquard, 
M. Steinhauser and C. Sturm for their collaboration and S. Uccirati for help in preparing this 
manuscript. 
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